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We re-calculate the exclusive semileptonic and nonleptonic decays of Be meson to a P-wave 
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^ I charmonium in terms of the improved Bethe-Salpeter (B-S) approach, which is developed recently. 
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^ . charged lepton spectrums for the semileptonic decays are precisely calculated. To test the concerned 
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^D . approach by comparing with experimental measurements when the experimental data are available, 
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'^'■ 

by some approaches else as well as the original B-S approach, which appeared in literature, are 
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I. INTRODUCTION 

The meson Be is the ground state of the double heavy (both of the components are heavy) 
quark-antiquark binding system (6c). In Stand Model (SM) it is an unique meson which 
carries two different heavy flavors explicitly, thus it decays weakly, that is very different 
from the ground states of flavor-hidden double heavy mesons, such as rjc and rji,. Namely 
Be decays via weak interaction (via virtual W emitting or annihilating) only, while the 
ground states of flavor-hidden double heavy mesons decay dominantly by annihilating to 
gluons (strong interaction) or/and photons (electronic interaction). The meson Be has very 
rich and experimentally accessible decay channels, so to study the decays of Be meson is 
specially interesting. By comparing the experimental and theoretical results of the decays 
of Be meson, we can also reach some insight into the binding effects of the heavy quark- 
antiquark system, which are of QCD nature, besides the knowledge of the weak interaction 
such as the CKM matrix elements etc. 

The meson Be was first experimentally discovered by the CDF collaboration at Fermilab 
through the semileptonic decay Be — ?► J/^j + I + u [i| , and soon it is confirmed not only by 
CDF itself via another decay channel Be -^ J /i}) + vr [2|, but also by the other collaboration 

n r 

DO at Fermilab [3]. The latest experimental report for its lifetime and mass in PDG |^ is 
Mb, = 6.277 ± 0.006 GeV and tb, = (0.453 ± 0.041) x lO'^^ g^ Because the cross section 
of Be production is comparatively smalL so to discover it is quite difficult in experiment. 
Whereas according to the estimates (5|-[7|, that LHC will produce about 5 x 10^°-Bc events 
per year, it is expected that more measurements of decays and production of the meson Be 
are available soon at LHC (LHCb, CMS, ATLAS), and it must push more studies of the 
decays of Be meson forward. So both experimental and theoretical studies on Be meson now 
become more interesting. 

In fact, the decays of Be meson can be divided into three categories: i). The anti-bottom 
quark b decays into c (or u) with c-quark being as a spectator; ii). The charm quark c decays 
into s (or d) with 6-quark being as a spectator; iii). The two components, b and c, annihilate 
weakly. According to the decay products we may realize which one or two even three of the 
categories play roles in a concerned decay, thus one can measure the CKM elements such as 
He, Vub, Vcs, Vcd through the decays. In the present paper, we are highlighting the decays 
of Be meson to a P-wave charmonium, and one may easily to realize that the decays being 
considered here belong to the category i). Since the lepton spectrum and the weak form 
factors, which relate to the binding effects (wave functions) precisely, may be measurable in 



semi-lepton decays as long as the experimental sample of decay events is great enough, so 
we will share quite a lot of lights on them. 

In fact, one may find a lot of theoretical methods to treat the semi-leptonic and non- 
eptonic decays of Be meson, such as the varieties of relativistic constituent quark models 



15J and QCD sum rules 
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ITJ I etc in the literature, and moreover one may realize 



that among the relativistic constituent quark models, the method presented in Ref. [18l| and 
adopted in Refs. J9l.ll0| is based on the instantaneous version 19| of the Bethe-Salpeter (B-S) 



equation 20|, and the 'instantaneous treatment' is also extended to the weak-current matrix 
elements using the Mandelstam formulation 2l|, while the adopted approach in Ref. |8| is 



different from the one presented in Ref. 18 



only in the kernel of the B-S equation and the 



'instantaneous treatment' etc. Recently in 22(| an improvement to that of 18| is proposed, 
and the relativistic effects in the binding systems and decays between the systems may be 
considered by the new development more properly, especially, considering the fact that, of 
the new development, the part (factor) for dealing with the binding effects has been applied 
to study (test) the spectra of positronium (a QED binding system) [23| and double heavy 



(see Refs. 



23 



flavor bin ding sy stems (QCD binding systems) 24 1 and quite satisfied results are obtained 

22| when experimental data are available 



2J| ) , so to test the new development 



in foreseeable future, in this paper we try to apply the development to the decays of B^ meson 
to a P-wave charmonium and to compare the obtained results with those obtained by old 
method in Ref. 18| and obtained by other theoretical approaches. Since we suspect that the 
decays of B^ meson to a P-wave charmonium might be more sensitive in testing the effects 
caused by the improvement than the decays of B^ meson to an S-wave charmonium, so here 
we focus our attention on the decays of B^ meson to a P-wave charmonium. 

The new development 22| contains two factors: one is about relativistic wave functions 
which describe bound states with definite quantum numbers, i.e. a relativistic form of wave 
functions (see Appendix C) which are solutions of the full Salpeter equation (see Appendix 
B). Note that here we solve the full equations Eqs. (B9, BIO, Bll), not only the first one 
Eq. (B9) as other authors did. The other factor of the improvement is about computing the 
weak-current matrix elements for the decays with the obtained relativistic wave functions 






n 



as input. It is more 'complete' than that as done in Refs. [9|, |lO|, ll8|, i.e. the 'complete' 
formula in Eq. (15). 

The paper is organized as follows: the formulations of the exclusive semi-leptonic and 
non-leptonic decays are outlined in Sec. II. The newly developed formulations, mainly for the 



matrix elements of the hadron weak decays, are presented in Sec. III. In Sec. IV, numerical 
calculations for the exclusive semi-leptonic decays and non-leptonic decays are described, 
the results and comparisons among the various approaches are presented. Finally the Sec. V 
is attributed to discussions. In Appendices, the formulations as necessary pieces for the 
calculations of the decays are given. 

II. THE FORMULATIONS FOR EXCLUSIVE SEMI-LEPTONIC DECAYS AND 
NON-LEPTONIC DECAYS 

Let us now derive the formulations for the exclusive semi-leptonic and non-leptonic decays 



precisely (mainly quoted from Ref. J22|) for numerical calculations later on. 

In the following subsections we will focus light on the matrix elements of weak currents, 
and show how to present the amplitudes of the semileptonic or nonleptonic decays via the 
matrix elements of weak currents precisely. In fact, one may see that the newly developed 
method mainly is about the matrix elements of weak currents. 

A. The semileptonic decays of Be meson 

The Fig. 1 is a typical Feymann diagram responsible for a semileptonic decay of Be meson 
to a charmonium. The corresponding amplitude for the decay can be written as: 
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FIG. 1: The Feynman diagram of a semileptonic decay of Be meson to a charmonium. 

where V^c is the CKM matrix element, {xc{hc){Pf)\J^\Bc{P)) is the hadronic weak-current 
matrix element responsible for the decay, and P, Py, py and pi are the momenta of initial 



state Be, the finial P-wave state of (cc) (i.e. h^ Xco, Xa, Xci and their excited states), the 
neutrino and the charged lepton respectively. 

Generally, the form factors are defined in terms of the matrix elements of weak current 
responsible for the decays appearing in Eq. ([1]). Namely for the decay of Be meson to scalar 
charmonium Xco, the form factors s+ and s_ are defined as follows: 

(x.o(P/)|^1i?c(P)) = 0, 

(Xco(P/)|^lSc(P)) = s^{P + P^Y + s4P- PfY . (2) 

For the decay of B^ meson to vector charmonium Xci-, the relevant form factors /, u\, U2 
and g are defined as follows: 

{Xci{Pf)\V^\Be{P)) = f{M + Mf)e^ + [u,P^ + u,P^f-j^, 

{Xci{Pf)\A^\Be{P)) = J^^;'"^'^'^PfPs^ ■ (3) 

For the decay of B^ meson to vector charmonium he, the relevant form factors Vq, V\, V2 
and V3 are defined as follows: 

{he{Pf)\V^\Be{P)) = Vo{M + Mf)e^ + [V,P^ + V^Pff-j^, 

{K{Pf)\A^\BeiP)) = Jl^^^^"""'^-^PpPf. ■ (4) 

For the decay of Be meson to tenser charmonium Xc2, the relevant form factors k, c\, C2 and 
h are defined as follows: 

(Xc2(P/)|Al5,(P)) = k{M + M^)e^°^ + ee.p^^{c,P^ + c^P^), 

1h P" 

{xAPsW-mp)) = m^m;-"*. 17 ^'""" w, . (5) 

In the case without considering polarization, we have the squared decay-amplitude with 
the polarizations in final states being summed: 
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^s.,s„S^.J.jT\'' = ^\V,e\\.h^^ , (6) 



where l^u is the leptonic tensor: 



and the hadronic tensor relating to the weak-current in Eq. ([T]) is 

h^-^ ^ j:s,^^^jB,iP)\j^+\xciK)iPf)){xciK)iPf)\r\B,iP)) 

= -ag'"' + /3++{P + PfYiP + PfY + /3+_(P + PfY{P - PfY 

+f3^+{p - PfY{P + PfY + /3_-(P - PfY{P - PfY 

+i^e^'-p-(P + p^)^(p - p^)^, (7) 

where the functions a, /3++, /3+-, /?-+, /3 , 7 are related to the form factors and we put 

the relations in Appendix A precisely. 

The total decay width F can be written as: 
1 f d^Pf d^pi d^p 



{27rY6\P-Pf-pi-p.)j:s.,s„s^,jT 



where Ef , Ei and E^ are the energies of the charmonium, the charged lepton and the neutrino 
respectively. If we define x = Ei/M, y = (P — PfY/M^, the differential width of the decay 
can be reduced to: 



+/3- 



++ 



d^T 
dxdy 



,,, ,oGlM^ [ 2a , 



647r3 M2 



M2- 



4(2x(l-^ + y) 



4x' 



+ (/5+_ + /?_+) 



M2 



2-4x 



V 






rrii 
M2 



27Z/ 



fi-^ 



M2 






(9) 



here M is the mass of the meson 5^, Mf is the mass of the charmonium in final state, and the 
total width of the decay is just an integration of the differential width i.e. T = j dx j dy-^^. 
Thus the key problem for calculating the semileptonic decays is turned to calculating the 
hadronic weak-current matrix elements. 

B. The nonleptonic decays of B^. meson 



In this subsection we mainly consider the nonleptonic two-body decays to a P-wave 
charmonium, i.e. decays Be — t- M1M2 where Mi is a P-wave charmonium and M2 is a 
common meson. Fig. 2 is the Feynman diagram for the decays via the relevant effective 



Hamiltonian H^ff 

Heff 



25 



26|: 



S^" 
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FIG. 2: The Feynman diagram of a nonleptonic decay of B^ to two mesons Mi (a charmonium) 
and M2 (a common meson). 

where Gp is the Fermi constant, q = d,s, Vij are the CKM matrix elements and Ci{fi) are the 
scale-dependent Wilson coefficients. Oi are the operators constructed by four quark fields 
and have J'^J^ structure as follows: 

of = [Vud{daUa)v-A + Vus{SaUa)v-A + Vcd{daCa)v-A + Vcs{SaCa)v-A]{Cfibfi)v-A , 
of = [Vud{daU^)v-A + Vus{SaUi3)v-A + Vcd{,daCp)v^A + Ks(SaC/3)y_A] (C/36a)y_A , 

O3 = (ga6a)y-A^(g^g^)y-A, O4 = (g/36a)y_A^(glg^)y-A, 

q' q' 

O5 = {qaba)v~Aj2(^Wl3)v+A, Oq = {qpba)v-Aj2(^'a(l'lB)v+A, 
q' q' 

3 3 

Or = ^iqaba)v-Aj2^<]'WlB)v+A, Og = -{qi3ba)v-AY. ^<l'(^Wv+A, 

<?' g' 

3 3 

O9 = ^{qaba)v-Aj2^Q'(^'l3l'l3)v-A, Oio = -(g/3fea)y-^ J] V(^a?/j) V-A, (H) 

q' q' 

where {qiq2)v-A = ^i7^(l — 75)q'2- The operators Oi and O2 are the current- current (tree) 
operators, 03,...,0q are the QCD-penguin operators and Oj,...,0io are the electroweak 
penguin operators. Since we calculate the decay up-to leading order, we just consider the 
contribution of Oi and O2. 

Here we apply the so-called naive factorization to H^jf i.e. the operators Oi [27|, so the 
nonleptonic two-body decay amplitude T can be reduced to a product of a transition matrix 
element of a weak current (Mi| J^|i?c) and an annihilation matrix element of another weak 
current {M2\J^\0): 

T = {M,M2\H,ff\B,) ^ {M,\J^\B,){M2\J,\0) , (12) 

while the annihilation matrix element is relating to a decay constant directly. The reason 
why we adopt the naive factorization here is that it works well enough due to the fact that 



all the decays concerned in this paper are 'constrained' to those in them the quark c as a 
'spectator' goes from initial Be meson into the final meson Mi always, thus as pointed by 



the authors of 28|, |29|] , in the concerned cases the corrections to the naive factorization are 
suppressed. 

Since Mi = Xc{hc), the matrix element {Mi\J^\Bc) is just the hadronic weak-current 
matrix element appearing in the previous subsection, but different from it by momentum 
transfer being fixed (owing to the decays are of one to two-body). The annihilation matrix 
element (M2I J^|0) with J^ = {qiq2)v-A is related to the decay constant of a 'common meson' 
M2 and can be measured via proper processes generally. 

Precisely, let us now 'restrict ourselves' to analyze the Be nonleptonic decays to the P- 
wave charmonium and the vr"*", p"*", etc, which are governed by the weak decay b — )■ cud, or to 
the P-wave charmonium and K~^, K*, etc, which are governed by the weak decay b — )■ cus. 
As an example, under naive factorization, we have the decay amplitude of B^ — )■ XcoP^ as 
follows: 

T{Be ^ xcop"-) = ^KfeV;>i(/i)(xco| Jl5c)(p+| J^IO) , (13) 

here ai = ci + -^€2 and N^ = 3 is the number of colors. 

Since (M2I J^|0) is relating to the decay constant of the meson M2 directly, so to calculate 
the widths of the non-leptonic decays is straightforward when the weak-current transition 
matrix elements {Mi\J'^\Bc{P)) are well calculated. Thus one may see that the problem 
to calculate the non-leptonic decays is essentially attributed to calculating the hadronic 
weak-current matrix elements (Mi|V^'^|i?c(-P)) and {Mi\A'^\Bc{P)) appearing in the above 
subsection for semileptonic decays. 

III. COMPUTATION OF THE TRANSITION-MATRIX ELEMENTS FOR 
WEAK-CURRENTS 

From the section above, we can see that to calculate the weak currents matrix elements 
{Mi\J^\Bc{P)) is the key problem for the concerned semileptonic and nonleptonic decays, 
so let us now explain the reason why and show how to apply the newly developed method 
|22| to calculate the matrix elements. In fact it is also to prepare necessary formulae for 
final numerical calculations. 

Here the weak-current matrix elements are for 'transitions' from a state of a double heavy 
meson to another double heavy meson. Due to the mass difference of the two states, the 
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relativistic effects for tlie transitions are great, tliat a proper formulation to deal with the 
relativistic effects is desired. It is known that the approach of relativistic B-S equation for the 
bound states and Mandelstam formulation for the transition matrix elements may be taken 
into account quite well, and furthermore the B-S equation and Mandelstam formulation 
even under 'instantaneous approximation' still works, because here the involved mesons are 



double heavy. While the newly developed method 22|, which applies the 'instantaneous 



approximation' to the current matrix elements and B-S equation completely, should be 



better than the original one in Ref. 18|, where the 'instantaneous approximation' is applied 
incompletely. The 'completeness' here means to apply it to the B-S equation, the solutions 
(B-S wave functions) and the transition matrix element (under Mandelstam formulation) 
properly, and let us outline it below. 

According to the Mandelstam formulation [21], the corresponding hadronic matrix ele- 
ments of weak current between the double heavy meson Be in initial state and the double 
heavy meson Xc{hc) in final state, appearing in Eq. ([1]), Eq. f fT2|) and Eq. f fT3|) . can be written 
as: 

(Xc(M(P/)|Jli?c(P)) 



-Tr 



X,,^^H.){P',q'){A-m,)x,AP,q)V,bY{l~lMPi-p'i) 



= tJ^Jr[x^^^^^^iP',q')ia, /'+ ^- mi)x,^(P, g)K.7^(l - 75)] , (14) 

where pi = aiP + q (ai = ^"V^ ), P2 = ^2-^ — Q («2 = rn^m ) ^^^ ^^^ momenta of c-quark 
and 6-quark respectively inside Be meson; p[ = a[Pf + q' {a[ = J^^m' )' P'2 — '^2-^/"?' ('^2 — 
rrP2m' ) ^^^ ^^ momenta of c-quark and c-quark respectively inside the P-wave charmonium 
Xc{hc)', moreover, for the final result (the last line of Eq. ( 1T^ ) we have P = Pf +pi +Py and 
q' = aiP + q — a'lPf. 

The newly developed method 22| essentially is to apply the 'instantaneous approxima- 
tion' to the current matrix elements and the B-S equation completely, to outline it and for 
'applying the instantaneous approximation' in a covariant way, we need to decompose the 
relative momentum q into two components: the time-like one glf and the space-like one qj_ 

as follows: 

P ■ a 

y ~ yy + y±' ^\\ — yi|2 ' yj- — y yji > 
(P ■ P') 

p'M -- p'M _|_ p'M p'M ^= 1 i pM p'A« ^ p'M P'/^ii • 

-T -T II -r -T _L , J- II _ ^2 -< ) -T ± — -T -* II ) 



and 



t'M 



fa I rti 

<1\\ + 9_L > 



la 
(1± 






g- = g,- + gl, g[f = (P • q'/M')P>^, 

where M is the mass of the meson B^., and we may further have two Lorentz invariant 
variables qp = ^ and qr = 



-ql- 



The 'instantaneous approximation' applying to the matrix element is just to carry out 
the integration of dqv. by a contour one on Eq. (IT^ precisely and to obtain the result below: 



{Xc{hc){Pf)\J'\B,{P))=i 

d^q± 

(27r)3 



d^q 



(2vr) 



Tr 



Xxc{m(^''^')(«i /'+ ^-mi)XsAP,(l)VcbY{l-l5] 



Tr< 



v'^H<i'±)^v^H<i±) + v^H<il)^r~iQ±) 



M^v^Hqi 



r 



where: 






^ +(gj 



^+-(g±) 



+</'' *^^^)m^ ^*^^^) 



A^(g±)r7(g±)A^(g±) 

M — cJi — a;2 

Ar(g±)^(9±)A2~(9±) 

M + cJi + u;2 

Ar(9p±)^(9±)A^(gpx) 

M - tJ2 - ^^2 - ^/ 

A+(gp±)r7(gx)A2~(gpx) 

M + W2 + CU2 - ^/ 



^'" 






75 J 



(15) 



-^'^n^pJ 



'+/^/ 



V 



^2^WpA.)nWpA.)K^{^_ 



PL) 



Ef 



UJ[ 



U2 



f' iQ'p±) 



K-iq'pL)viq'p±)K-{q', 



p±) 



i,'-+{q' 



P±) 



^'^~iq'p±) 



Ef + uj[+ uj'2 

K{q'pLW{q'p^)A'^{q'p^) 

M -UJ2-0J2- Ef 

K'-{q'p±)v'{q'p±)K~{q'pL) 



(16) 



M + UJ2 + UJ2- Ef 

V'*''(Q'±)) i^^Kq^) ^^d (^'*-'(gp_|_), '?/''*"' (gp_|_) are B-S wave functions as the B-S equation solu- 
tions under 'complete instantaneous approximation' 23|, [2^ and with 'energy projection' 
A^ of the mesons in initial and finial states properly. The precise definitions of the 'energy 
projection' and the B-S 'vertex' rjp, fjp {rj' p, rj'p) are presented in Appendix B. One may 
also see that the four equations, Eqs. flB9[IB10llBlip . are B-S equations under the complete 
instantaneous approximation, instead of the incomplete instantaneous approximation which 
only considering the Eq. ( lB9p . 

Namely the 'improvements' from the 'newly development method' are attributed to: i). 
with the complete instantaneous approximation to current matrix element, as a result, there 
are six terms in the squared bracket of Eq. (llSp instead of the first term 



(Xc(M(P/)|Jl5c(P)) 



|^Tr{^--(gl)^^+-^(g,)7'^(l-75)} 



(17) 
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is only kept; ii). the B-S wave functions hidden in ip^^{q±), 4'^^{q±) and <^'*-'(q'pj_), 4''^-' {(l'p±) 
are solved under complete instantaneous approximation to the B-S equation. For the point 
i), since the considered double heavy meson, Be or Xc{hc), is weak binding system i.e. the 
binding energy e = M — Ui — uj2 (or £ = Ef — uj[ — ijJ2) is small {jj <^ 0(1)), thus 
from Eq. (fT6l) we are sure that if~^~^{q±) and 'f'^^iq'p±) are much greater than the others 
(p^^{q_i_), ■ip'^^{q±) and Lp'^^{q'pj_), 4'''^^{q'p_i), so that using the Eq. fITTj) instead of Eq. flTSj) is 
a very good approximation, which we have precisely examined by considering the decay 
Be -^ Xco^^i as an example: in fact, the contributions of the second term and third term of 
Eq. (TT5|) to the form factor are less than the one of first term of Eq. (fTSjl roughly by a factor 
10^^ ~ 10^^ times. If the first three terms are considered, the decay width is 1.85 x 10~^^ 
GeV, while if only the first term is considered, the decay width is 1.87 x 10^^^ GeV, i.e. the 
two results are very similar. So the approximation is very good and we may use Eq. flT7|) 
instead of Eq. (fT5|) to compute the weak-current matrix elements safely. 

IV. NUMERICAL CALCULATIONS AND RESULTS WITH PROPER COMPAR- 
ISONS 

In this section, based on the formulations obtained in the paper, we evaluate the decay 
widths for semileptonic and nonleptonic decays and some interesting quantities else for 
semileptonic decays, such as form factors and charged lepton spectrum etc and then discuss 
them briefly. 

First of all, we need to fix the parameters appearing in the framework. We adjusted 
the parameters a = e = 2.7183, A = 0.21 GeV^, Aqcd = 0.27 GeV, mt = 4.96 GeV, 



24 



30, 



3l| . which as the best 



rric = 1.62 GeV and Vq for the B-S kernel as those in Refs. 
input for spectroscopy, then the spectra of the mesons and the masses Mp^ = 6.276 GeV, 
M^^o = 3.414 GeV, M^^, = 3.510 GeV, M^^^ = 3.555 GeV, M^^ = 3.526 GeV etc ^2^], which 
are used in this paper, are obtained, moreover the decay constants, average energies as well 



as annihilations of quarkonia are fitted 30l-l32|. 



With the obtained B-S wave functions (under the formulation defined in Appendix B) 
and as a next step, we substitute the functions into '^~^^{qi_) and i^''^~^{q'pj_), so that they are 
related to the components of the B-S wave functions precisely as depicted in Appendix C. 
With the formula Eq. flTTl) . finally we represent the hadronic transition weak-current matrix 
elements as proper integrations of the components of the B.-S. wave functions. As final re- 
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TABLE I: The semileptonic decay widths (in the unit 10 GeV) 



Mode 


This work 


[12j 


[13] 


[15] 


[10] 


[16] 


[17] 


Bt -^ Xcoei^ 


1.87 ±0.46 


1.27 


2.52 


1.55 


1.69 


2.60±0.73 




Bt -^ XcQ-riy 


0.23 ±0.12 


0.11 


0.26 


0.19 


0.25 


0.7±0.23 




Bt -^ Xciev 


1.52 ±0.45 


1.18 


1.40 


0.94 


2.21 


2.09±0.60 




Bt -^ Xciry 


0.14 ±0.10 


0.13 


0.17 


0.10 


0.35 


0.21±0.06 




Be ~^ Xc2ev 


1.50 ±0.39 


2.27 


2.92 


1.89 


2.73 






Bt -^ Xc2TV 


0.12 ±0.07 


0.13 


0.20 


0.13 


0.42 






-B+ — )• hceiy 


3.98 ±1.10 


1.38 


4.42 


2.4 


2.51 


2.03±0.57 


4.2±2.1 


B^ — > hcTv 


0.28 ± 0.20 


0.11 


0.38 


0.21 


0.36 


0.20±0.05 


0.53±0.26 



suits of this paper, the decay widths for the semileptonic and nonleptonic decays and some 
interesting quantities else for the semileptonic decays, such as form factors and charged lep- 
ton spectrum etc, are straightforwardly calculated numerically. In the following subsections 
we present the results for the semileptonic decays and nonleptonic decays separately. 

A. The semi-leptonic decays 

When the weak current transition matrix element for a definite semi-leptonic decay is 
calculated precisely and the values of the CKM matrix elements \Vud\ = 0.974, \Vus\ = 0.225, 
I Vfecl = 0.0406 j^ are given, not only the decay width can be calculated straightforwardly, but 
also the form factors may be extracted out. Moreover as 'semifinished product', the spectrum 
of the charged lepton which may be measurable experimentally can be also acquired too. 

Namely the functions a, /3++, /3+_, /3_+, /3 , 7 appearing in the spectrum of the charged 

lepton (see Eq. (Q) are related to the form factors directly as shown in Appendix A precisely. 
Therefore when we calculate and present the results for semi-leptonic decays, not only those 
of the decay widths but also the spectrums of the charged lepton in the decays are considered. 
Since r lepton is quite massive and m^ ~ rrie is quite a good approximation for the Be meson 
decays, so when we calculate and present the widths and the spectrums of the charged lepton 
for the decays, only the cases that the lepton being electron or r are considered. 

Note that since the input B-S wave functions by solving the B-S equation for the double 
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FIG. 3: The form factors of the B^ decays to a P-wave charmonium defined as in Eq. ljAip . Eq. ()A2p . 
Eq.dMl) and Eq.dAH) and t = q^ = {P - Pjf = M'^ + M'j - 2MEf {tm is the maximum of t). 



heavy mesons which are involved in the transition matrix elements of weak current have 
uncertainties, due to the parameters fitting to fix the B-S kernel and quark masses, the way 
to solve the B-S equation numerically, and the approximation from Eq. flTSj) to Eq. f lT7|) 
for the transition matrix elements of the weak currents is taken etc, so in the numerical 
results obtained finally there are certain errors. To consider the uncertainties caused by 
the input parameters, we changed all the input parameters simultaneously within 5% of 
the center values, then we get the uncertainties of numerical results for the semi-leptonic 
decays and the non-leptonic decays shown in Table. I. We find that the uncertainties of the 
decays Be — ?■ hdxc) + e + Ue vary up to 30% of center values, while the uncertainties of 
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FIG. 4: The energy spectrums of the charged lepton in the Be semileptonic decays to P-wave 
charmoniums. The left figure is for B^ — s- Xc0,i,2(^c)ei^ and the right figure is for Be — s- XcQ,i,2{hc)'''v. 
Where the sohd hues are the results for XcO, the dash lines are for Xci^ the dot lines are for Xc2 and 
the dot-dash lines are for he- 



Be — )■ hc{xc) + T + Ut are up to 60% in Table. I, the reason is that the phase spaces for 
Be — 7- hcixc) + T + Vr are smaller than the ones for Be — ?■ hdxc) + e + Ue because of the 
heavy r lepton, and the the uncertainties for the former are more sensitive to the changes 
of the phase space than the latter. 

To compare with the results obtained by the other approaches, we present the decay 
widths calculated out this work with error bar and the results obtained by the other ap- 
proaches by putting them together in a table i.e. Table I. 

In addition we also present the obtained form factors and the spectrums of the charged 
lepton in the decays in Fig. 3 and Fig. 4 respectively. To compare with the results of the 
previous work Ref . [10] , we draw the curves of the spectrums of charged lepton obtained by 
this work and the work Ref. 10| in Fig. 5. Whereas in order to see the tendency of the form 
factors and the lepton spectrum clearly and we suspect that at present stage it is enough, so 
in the figures we draw the curves with the center values but not involve the errors precisely. 
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FIG. 5: The energy spectruins of the charged lepton in the Be semileptonic decays to P-wave 
charmoniums respectively. The sohd hnes are the results of this work, the dash lines are the results 
of 
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B. The non-leptonic decays 



The exclusive non-leptonic decays are of two-body in final states, thus the hadronic 
transition matrix elements of weak-currents appearing in Eq. (fT3|) have a fixed momentum 



transfer t 



m 



M2 



(the mass squared of the other meson M2 in the decay Be — t- M1M2 and 



Ml = Xc or he). In fact the transition matrix elements have been already calculated in the 
above subsection of semi-leptonic decays. To calculate the decay widths, from Eq. flT3|) . now 
we need to calculate the annihilation matrix element of the weak current such as (M2I Jju|0) 
additionally. It is known that the annihilation matrix element is related to the 'decay 
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TABLE II: The nonleptonic decay widths (in the unit 10 GeV) 



Mode 


This work 


[12] 


m 


m 


m 


Bt -^ XcO^+ 
Bt ^ Xcivr+ 
Bt ^ Xc2vr+ 
Bt -^ K-K+ 


(0.34±0.04)af 

(0.023 ± 0.002)af 

(0.24±0.05)a? 

(1.10±0.16)af 


0.23af 
0.22af 
0.41af 
O.Slaf 


0.622af 

0.076af 

0.518af 

1.24af 


0.28af 

O.OOlSaf 

0.24af 

0.58af 


0.317a? 
0.0815a? 
0.277a? 
0.569a? 


Bt -^ XcOP+ 
Bt ^ XciP+ 
Bt -^ Xc2P+ 
Bt -^ hcp-^ 


(0.85±0.10)af 
(0.25±0.02)af 
(0.62±0.12)af 
(2.50 ± 0.50)af 


0.64af 
0.16af 
1.18af 
l.llaf 


1.47af 

0.326af 

1.33af 

2.78af 


0.73af 
O.llaf 
0.71af 
1.41af 


0.806a? 

0.331a? 

0.579a? 

1.40a? 


Bt ^ XcoK+ 
Bt ^ XciK+ 

Bt ^ Xc2i^+ 

Bt -^ hcK+ 


(0.026 ± 0.003)af 

(0.0018 ±0.0002)af 

(0.018 ± 0.003)af 

(0.082 ±0.012)af 


0.018af 
O.OlGaf 
0.031af 
0.039af 


0.0472af 
0.0057af 
0.0384af 
0.0939af 


0.022af 

0.00012a? 

0.018af 

0.045af 


0.00235a? 
0.0058a? 
0.00199a? 
0.0043a? 


Bt ^ XcoK*+ 
Bt ^ Xcii^*+ 

5+ ^ Xc2i^*+ 


(0.050 ± 0.006)af 

(0.018 ± 0.001)a? 

(0.037 ± 0.007)af 

(0.14±0.02)af 


0.045af 
O.Ola? 
0.082af 
0.077af 


0.0787af 
0.0201af 
0.0732af 
0.146af 


0.041af 
0.008af 
0.041af 
0.078af 


0.00443a? 
0.00205a? 
0.00348a? 
0.0076a? 



constant' /mj directly, and the decay constant fp, fy or /^ of a pseudoscalar meson, a 
vector meson or an axial vector meson may be extracted from experimental data for the 
pure leptonic decays of the relevant mesons, but they may also be calculated by models, 
such as the one in Ref. |30| although there are some debates. In this work we adopt the values 
of the decay constants: /^ = 0.130 GeV, fp = 0.205 GeV, fx = 0.156 GeV, fx* = 0.217 
GeV etc for numerical calculations. Then the relevant decay widths for the concerned non- 
leptonic decays are calculated. As the final results, we present the decay widths by our 
method and the others' methods else in Table II. Note that the uncertainties in Table II are 
estimated as done in the previous subsection for semileptonic decays. 

For comparison precisely with the other approaches and experimental measurements in 
future, we take the values ai = 1.14 for non- leptonic decays as done in most references, and 
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TABLE III: Branching ratios (in %) of Be decays calculated for the Be lifetime tb^ = 0.453 ps and 
ai = 1.14. 



Decay 


Br 


Decay 


Br 


Bt -^ Xc^ev 


0.13 ±0.03 


Bt -^ Xcqtu 


0.016 ±0.008 


Bt -^ Xcieiy 


0.11 ±0.03 


Bt -^ XclTV 


0.0097 ± 0.0065 


Bt -^ Xc2eiy 


0.10 ±0.03 


Bt -^ Xc2Tiy 


0.0082 ± 0.0048 


B^ -> hccu 


0.28 ± 0.08 


Bt -> hcTU 


0.019 ±0.013 


Bt -^ Xcovr+ 


0.031 ± 0.004 


Bt ^ XcoP+ 


0.076 ± 0.009 


B+ -^ Xci^+ 


0.0021 ± 0.0002 


Bt -^ XciP+ 


0.023 ± 0.002 


B+ -^ Xc27r+ 


0.021 ± 0.005 


Bt -^ Xc2P+ 


0.056 ±0.011 


B+ -^ hcTT+ 


0.098 ± 0.015 


Bt -^ hcp"^ 


0.22 ±0.04 


Bt -^ XcoK+ 


0.0023 ± 0.0003 


Bt ^ XcoK*+ 


0.0045 ± 0.0006 


Bt ^ XciK+ 


0.00016 ± 0.00002 


Bt ^ XciK*+ 


0.0017 ±0.0001 


Bt ^ Xc2K+ 


0.0016 ± 0.0003 


Bt ^ Xc2K*+ 


0.0033 ± 0.0006 


Bt -^ hcK+ 


0.0074 ±0.0011 


Bt -^ KK*+ 


0.013 ±0.002 



the experimental value of Be lifetime tb^ = 0.453 ps as well, we calculate branching ratios 
of the decays and put them in Table III. 

V. DISCUSSIONS AND CONCLUSIONS 



In Sec. IV, the form factors (Fig. 3), energy spectrums of the charge leptons (Fig. 4 and 
Fig. 5), decay widths (Table I) for the semileptonic decays, and the decay widths for non- 
leptonic decays (Table II) are presented. Specially in tables some comparisons with other 
approaches else are also given. Thus one may read off a lot of interesting matters already. 

Since the form factors for the semi-leptonic decays, which are directly related to overlap- 
ping integrations of the components of the B-S wave functions of the initial and final states 
as shown in Appendix C, are comparatively difficult to be measured, so in Fig. 3 we show the 
behaviors of the form factors briefiy (without errors). Whereas the energy spectrums of the 
charged lepton in the decays may be measured not so difficult, as long as the event example 
is great enough and the abihties of the detector are strong enough, and to see the differences 
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between the spectrums of electron and r lepton clearly in Fig. 4 we plot the curves with 
center values without theoretical uncertainties. Moreover to see the differences between this 
work and the ones 



10| , in Fig. 5 we plot the spectrums of electron obtained by this work vs 
the ones [10| obtained by previous approach and for both of them only center value without 
theoretical uncertainties are taken. Since the spectrums of muon (/i) is very similar to that 
of electron in exclusive semi-leptonic decays, thus we do not present the spectrums of muon 
at all. From Fig. 4 we can see the difference in the energy spectrums among the Be decays to 
different P-wave charmonia clearly, although the results of electron is greater than the one 
of T lepton. From Fig. 5 we can see that the difference in the energy spectrums of electron 
due to different approaches: the difference caused by newly improved approach and by the 
previous approach can be quite sizable and can be tested experimentally in future. For 
the widths of the decays, from Table I and Table II, both the semi-leptonic decays and the 
non-leptonic decays, one may see that in general the results of this work fall into the region 
of the predictions by various models, but the distribution of the predictions is quite wide, 
so future experimental data will be critical and may conclude which one of the predictions 
is more reliable. 

Considering the fact that the substantial tests of the i?c-meson decays have not been 
started yet, although the meson B^ has been observed at Tevatron for years and LHC 
is running now, according to the estimates of the production at LHC, one may believe 
reasonably that the tests of the predictions on the Be decays will be started with LHC more 
measurements available. From theoretical point of view, we think that the newly improved 
approach works better than the previous one, this trust need to be tested by experiments. 



We would also like to note here that according to the estimates |33H36| of the production 
at an e^e+ collider running at CM energy \/S ~ mz [mz is Z-boson mass) with very high 
luminosity {L = 10^^~^^cm^^s~^) i.e. a "Super- Z- Factory" and considering the advantages, 
may be more suitable to test the approaches by measuring the decays precisely than that to 
do them at hadronic collider such as Tevatron or LHC, because at such a Super- Z-Factory 
numerous Be mesons may be produced and the energy- momentum of the produced Be meson, 
as the e~e~^ one of the collider, is precisely known in an e"*" — e~ collider environment. 
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Appendix A: The functions a, /3++, /3h , /? i-, /3 , 7 

Here according to the P-wave charmonium appearing in the final state we present the 

useful functions «, /?++, /?+-, /? h, /? , 7 how precisely to relate to the form factors in turn. 

a). When Be decays to Xco'- 

Since the matrix elements of weak currents are described in terms of two form factors 

(s+,s-): 

(Xco(P/)|Vlfic(P)) = 0, 
{Xco{Pf)\A^^\B,{P)) = s+{P + PfY + s^{P- PfY , 

then the functions are read as 

/3++ = 4, /3+- = /?_+ = s+s_, /3__ = s^ (Al) 

b). When B^ decays to Xci- 

Since the matrix elements of weak currents can be described in terms of four form factors 

(/,Ml,M2,^): 

,e-P 



{Xci{Pf)\V^\B,{P)) = f{M + Mf)e>' + [mP^^ + U2P 



f^ M ' 



{XcAPf)\A^\Bc{P)) = -M^;'''""'^'^PpPf.^ 



then the functions are read as 

f2 I aT\/T2^,2^ 



a = f^ + AM^gl0}. 



Jl_ /i#2„2„. , -L 



^-- = Tm-'"'^'^2 



2 












/l«+ + -7J2^U^ 



^2 1 1 /\/f z? M v 

T = -2/ig, (A2) 
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when setting f, = f{M + Mj),u+ = ^^^, u. = ^-^^, g. 



M+Mf' 



c). When Bf. decays to he- 

Since the matrix elements of weak currents can be described in terms of four invariant 
form factors (Vq, Vi, V2, V3): 



{K{Pf)\V^\B,{P)) = Vo{M + Mf)e^ + [ViP'^ + ^2^, 



^J M ' 



{h,{Pf)\A'^\B,{P)) 



2V. 



M + M 



-te^^'-e,P,Pf^ 



f 



then the functions are read as 



a = ft + AM-'gipj 



P- 



++ AM} 



^' M'gh + \ 



My ^' 



fia+ + 



Mj +' 



P+. = f3.+ = gfiM' - M') - 



M2^ 



ft I,., , IMEf^^ , ivrp". 

/i(a+ + a_) — -— — j-Ji(a,+ — 0,-) + a+a 



^' AM] 2 



2 M| 



M} ' 



,2, , fl (MEf , \ ^^ , 2MV/ 



/5- = ~9i{M' + 2ME^ + Mj) + ^- [-j^ + Ij /la- + a^ 



M| ' 



7 = -2/1^1 , 



(A3) 



when setting /i = \/o(^ + ^/), ^h 



2M : 



a_ 



{Vi-y2) 



VL 



2M ' 5*1 M+Mf • 



d). When Be decays to Xc2'- 

Since the matrix elements of weak currents can be described in terms of four form factors 
{k,ci,C2,h): 



{Xc2{Pf)\A^\B,iP)) = k{M + M;)e"^^ + e^^^{c,P^ + c,P^) , 



M 



M2 



(Xc2(P/)|V^1fic(P)) 



9h P" 



P^fa^ 



M + Mf ^M 
where ^^/^(e"^) is the polarization tensor of tensor meson, then the functions are read as 



2^2\ 



+-ckic+ 



a = -{kl + AMhipf 

(3++ = ^-^M'y+-c'cl 
'^^ 8M| 2 ^ 3 + 

'M\l-y) + MJ 1\ kl (M\l-y) + M] 1 



AM] 



AM] 



/3+- = /3-+ 



CrC-i C/t-| , o 



8M] ■ 2 



+ ^(M^-Mj) + 



6 



' M\l-y) + M] Y 1 
4M2 j 4 



-c^c+c_ 
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2 (M\l-y) + M] 1\ 2 (M\l-y) + Mj V 



^-f(2(M^ + M|)-M^,) + | 



+r^^^-i 4M^ 2J+y( 4M| +2J ' 

7 = -c/iifci , (A4) 

when setting c = ^, fci = /c(l + ^), c+ = ^, c_ = ^^, h = MiM+M^y 



Appendix B: The B-S equation under 'complete instantaneous approximation' 

In this appendix we outhne the 'complete instantaneous approximation' onto the Bethe- 
Salpeter equation when it has an instantaneous kernel, which describes a double heavy meson 
quite well. 

The Bethe-Salpeter equation |20| is read as 

{A - m^)Xp{q){^2 + m2)=i j -n^y{P. k, q)Xpik), (Bl) 

where Xpi^l) is B-S wave function of the relevant bound state, P is the four momentum of 
the meson state and pi, p2, rrii, m2 are the momenta and constituent masses of the quark 
and anti-quark respectively. From the definition, they relate to the total momentum P and 
relative momentum q as follows: 

7-. "^1 

Pi = aiP + q, «i — 



P2 = a2-P - q, Oi2 



mi + 7712' 
m2 



rrii + m2 

The interaction kernel V{P, k, q) for a double heavy system, being instantaneous approx- 
imately, can be treated as a potential after d oing instantaneous approximation, i.e. the 



kernel take the simple form (in the rest frame) 19 1 



V{P,k,q)^V{\k~(^) 



For various usages, we divide the relative momentum q into two parts. 



q' = q^ + ql, < = 1^^"' ll^l'-^l 
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where M is the mass of the meson, and we may have two Lorentz invariant variables: 

For the convenience below, let us introduce the definitions: 

^M'D^^I^XA^I.^'I). v{q'^)^l^,V{k^,q±)Mkl), (B2) 

then the B-S equation can be rewritten as 

X(9||, g±) = Si{pi)r]{q±)S2{p2) ■ (B3) 

Owing to Eqs. flB2t IB3p . it is reasonable and for convenience we may call i]{q±) as 'instan- 
taneous B-S vertex'. The propagator of quark or anti-quark may be decomposed: 

A+(g±) , Ar(g±) 



Si{pi 



J{i)qp + ttjM — Wi + it J{i)qp + ajM — Wi + ie 



where i=l, 2 for quark and anti-quark respectively, and J{i) = (—1)*+^, ui = Jml + g|,. 



UJ2 = \/Tn2 + qx, and A^ , Ag are the generalized energy projection operators. 



and have the properties: 



^Uq'pJ + Kp^^pJ = T-, ' A±,(g^J^Afp(g^J = 



A±p(g^J^A±p(g^J = A=^p(g^J , (B5) 

The instantaneous approximation to the B-S equation is to do contour integration over 
qp on both sides of Eq. fIBSp . and obtains: 

^ (^^) ^ A^(g±)r7(g^)A+(g^) _ Ar(g±)r/(gx)A2-(g^) ^^^^ 

If we introduce the notations: 

^J±(g^) ^ Af (g^)^^,(g^)^A±(g^) , (B7) 

we have 

<^p(g±) = <+(g±) + <"(g±) + ^;+(g±) + ^;"(g±) , (B8) 
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With the properties Eq. ( ]B5[) and notations Eq. (lB7p . the full Salpeter equation Eq. f lB6P 
can be written as 

(M-u,- u;2X+(g±) = Atiq±)v{q±)^tiQ±) , (B9) 

{M + UI+ U2)if;-{q±) = -A^{q±)r]{q^)A^{q±) , (BIO) 

^;-iq±)=^;'-iq±) = 0. (Bll) 

The normalization condition for the B-S equations now is read as: 

The couple equations Eq. (IB9p . Eq. flBlOp and Eq. fIBlip with the normalization condition 
Eq. (1B12P are the final B-S (Salpeter) equation under 'complete instantaneous approxima- 



tion' vs the previous one i.e. Salpeter equation [19| where only Eq. ( JBOP is considered. 

In addition, note that in the model used here for the double heavy quark-antiquark 
systems, the QCD-inspired interaction kernel V, being instantaneous approximately and 
dictating the Cornell potential which is composed by a linear scalar interaction plus a vector 
interaction, is read as: 

Vs{q) = -{- + Vo)6%q) ^ ^ 



a Tv"^ {(p + q;2)2' 

as{q) 



yM = -—,7^I^v (B13) 



where the QCD running coupling constant as{q) = 33_2jv lo (a+^/A^ — )' ^^^ constants 
A, a, a, Vq and Aqcd are the parameters characterizing the potential. 



Appendix C: The reduced wave functions (p~^~^{q) and the form factors 

In the appendix we present the reduced wave functions '^'^'^{q) (and 'ip~^~{q}) which 
directly relate to the solutions by newly solving the obtained coupled equations Eq. (lB9p . 
Eq. flBlOp and Eq. fIBlip under a new approach. The key point of the new approach is to 



solve the B-S equation according to the quantum numbers of the concerned bound states 



respectively [24 



30 



3l|, i.e. to solve the equation under the new approach we need to give 
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the most general formulation for the wave function first. Therefore for the present usage, 
in this appendix, we precisely quote the solutions for the low-laying bound states Be meson 
with quantum numbers J^ = 0^, Xco with quantum numbers J^*" = O"*""^, Xa with quantum 



numbers J^*" 
J'^'^ = 1+" from 



I++ 



24 



y£2 
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with quantum numbers J^'-' = 2"^+ and he with quantum numbers 



31| . and then we write down the reduced wave functions '^^^{q) 



and the form factors accordingly. 

When the weak-current matrix elements are computed precisely, as an intermediate step, 
the form factors can be represented as overlapping integrations of the components appearing 
in the B-S solutions, thus in this appendix we also give the formulas of the form factors in 
terms of the 'overlapping integrations'. 
a). For Be meson with quantum numbers J^ = 0~ 

The B-S wave function (solution of Eq. flB9p . Eq. flBlOp and Eq. fIBlip of B^ meson with 
J^ = 0~ is read as: 



^bM = m 



M 



/i(9l{l- 



^±{wi + W2) 
m2Wi + miW2 -' 



\ 



+ /2(g){i + 



4±{W2 -Wi) ■ 

miW2 + ni2Wi 



75, 



(CI) 



where M, P are the mass and the total momentum of the meson Be-, q± = (0,g), g* is the 
relative momentum of quark and anti-quark in the meson, so qj_ = —q"^. 
Then we can rewrite the reduced wave function: 



^t:{^ = bi 



Be 



^2 + — + &3 



+ h 



M 



75, 



(C2) 



where 



hi 



M ( mi + rn2 

2 \ Wi + W2 

(mi - 1712) 



•miW2 + m2Wi 
In Appendix. B in Eq. (lB2p . we have 



, W1+W2 

62 = ■ : 

mi + 1712 
(wi + W2) 
{miW2 + m2Wi) 



v{Qa 



d^kV{k)M 



M 



/i(^){l 



_l(wii +W21) I 



+AW{i+ fa''""-"'"' } 



miW2i + m2Wii 



7712^11 + miW2l 

75, 



(C3) 



where wu = J ml — k\, W21 
According to Eq. (JCll) . 



■ni.(ii 



ml - kl, V{k) = Vsik) + K(A;)7° ® 7" 



rf3A;(K(fc) + K(fc)7°®7° 
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M 



M '- m2Wii + miW2i ' ^ miW2i + m2Wii 



75 



M 






75- 



(C4) 



91 = j d%Vs - K]/i(^), 92 = j d?k[Vs - K]/2(^), 



\q\^ miW2i + m2Wii J 

So we can also write down the wave function of ip^ {<l±)i 



(Wii + W2l) 

miW2i + m2Wii 



ip^ (gj 



M + UJ2 + UJ2-Ef 






75- 



(C5) 



set tt 



4:WlW2{M+LLJ2+i02—Ef) 



where the symbol ', denotes of the final state, and 



ni = tt[giM{-q^ + mim2 - ^1^2) + 92M{m2Wi - miW2) + 9ziwi + ^2)1?^ + 9iij^i + "^2)1?^], 

^2 = tt[giM{m2Wi - miW2) + 92M{q'^ + mim2 - W1W2) + 9zij^i - m2)q^ + 9i{wi - 1^2)?^], 
^3 = tt[-giM{wi + ^2) - 92M{mi - 7112) + gz{q^ + mim2 + W1W2) + gi{m2Wi + ■miW2)], 
Hi = tt[giM{mi + 1112) + g2M{wi - W2) - g-i{ni2Wi + miW2) - 9i{-<f + ^17712 + ^1^2)]- 

b). For the charmonium XcQ {J^'" = 0+^) and the form factors s+ and s 

The B-S wave function (solution of Eq. ( JB9I) . Eq. ( JBIOI) and Eq. ( JBllI) under new method 
to solve the coupled equations) of Xco is read as: 



^xJq') = fiw) 41 + f2W)^ + ami + f'M) h 



Mf 



(C6) 



with constraints on the components of wave function, for the charmonium, m'^ = m'2,w[ 
w'2, we get: 



m) 



fi{q') = 0. 



Mfm[ ' 

where Mf, Pf are the mass and the total momentum of final meson Xco, ?! = (0, q'), q' is the 

_, 2 
relative momentum of quark and anti-quark in the meson, so q'1 = —q' . Then the reduced 



wave function '^tp{q') as: 



Vxli^') = «i 



^l + .2^ + a3-.a4 



(C7) 
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with 



«.4 /;tf) + «?)^ 



w\ 



02 



TTl 



I ' 



as = -^, 04 = 0. 



m 



The wave function of ip' ^(q'px) is 



K-Wp±W{q'p^)K'-{q' 



p±) 



M -UJ2-UJ2- Ef 



n. 



no , , +% + n4 



M 



/ 



w '« 



Set tt' 



4M)5^^(M-tJ2-^--B/) 



where 



and 



n\ = tt'[-2gW^ + 2g':,Mfm[], n'^ = 0, 
^3 = tt'[-2g[m\q''^ + 2g'^Mfm% n'^ = tt' [-2g[w[q''^ + 2g'^M fm'^w[] 



g, = I d'k'iV, - K]^/((fc'), g'2 = fd'k'lV, - V^jtj^f^ik'), 
\q \ \q \ 



g', = Jd^k'[Vs + V, 



f[ik')k " 
Mfm[ 



^ g', = 0. 



I ' 



With Eq. ( TT71) . the form factors may be presented by overlapping integrations: 
1 f d^q Aaibi 



2 J i2TifMMf 
,2 



azh2Mf + auEf{a2b2Ef + Mf + a2hq- Pj] 
+b3{Mfq^ + aiiMfq- Pf) + M{a2hq^ - anas^a^/ - aiiMj) 

+M- 



^qcos _^^ _ _j_^^^^^^^^ _ ^^^^j^^ + Mf + 0264? ■ Pf) 



IP, 



/I 



M- 



(C9) 



1 f d^q 4aibi 



2 J (27r)3 MM/ 



asb2Mf + aiiEf{a2b2Ef + Mf + a2b^q- Pf 



+b3{Mfq'^ + aiiMfq- Pf) - M{a2biq^ - aiia2&2-£^/ - aiiM/ 

-M^^^il + -^){a2b2Ef - asb^Mf + Mf + 0264?- Pf) 

\Pf\ M 



rrij^+mj 



(CIO) 



where an = a[ 

c). For the charmonium Xci {J^^ = 1^^) and form factors the /, Ui, U2, g 

The B-S wave function (solution of Eq. (IB9p . Eq. (IBlOp and Eq. (IBlip under new method 
to solve the coupled equations) of Xci is read as: 

^UQ')=^^,^c.^Pf^Te^[f'ii^')MfY+f2iq') ^fY+fsiQ') i'^r+^aq')e^'"''Pf.q'^,l8l,/Mf]/M'f 

(Cll) 
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where e is the polarization vector of axial vector meson and with the constraint on the 
components: 



m\ 



Then the reduced wave function iptp^{q') as: 



2 



(C12) 



with 



2 \ nil J "^1 ""^1 



With Eq. ( !T7|1 . the form factors may be presented by overlapping integrations: 



/ 



d^q 



Aaibi 



{27iy MVM + Mf 



(ciiM'jq^ - a2biM'jq^ - a2b2Efq- Pf 



+ (a4 - a2b4)iq- Pff + al^a^EJPf + a2b2aii{E) - %Mj) - 2ana4EJq- Pf 



+ana2b4EJq- Pf + b^EfMf{q' - anq- Pf) ) + y (cos' ^ - 1) (M;(a4 - 0364) 
-\-b2.Ef{Mf + a^q- Pf - aua^Pf 



(C13) 



Ml 



(Pq AaibiM 
(27r)3 M| 



aiiE 



f 



M2 



(a2fe2Mj + Mfbsq- Pf 



+a4{anEfM] + foslMjg^ + (g-. p^)2 _ c^^.g-. PfEJ))) - «n^J-^''°'^ 



M2|P^ 



/I 



'04 — 02^4) M 



/ 



+63£'/(— aiia4-E? + Mj + 04^- P/ + aua^Mf 



Ef qcos9 



f\ 



a2b2MJ + Mfb^^q- Pf +ai{aiiEfMJ + b^{MJq' + (g ■ PfY - a^q- PfEJ)) 



q' 



2M2|Pf|2 



/I 



-M^ + (2E; + M;) cos' ^) f M;(a4 - 0264) 



+b'iEf{Mf + a4g ■ Pf - aua^Pf 



(C14) 



«2 



d^q AaibiM 
(27r)3 Ml 
EfqcosO 

^ M\Pf\ 



— {b^Mfq^ + (a2&2 + «iia4-E'/)(aii£'^ - g ■ Pf 



M 



auMfbsEf + (a2&4 - 04) ((f- -P/ - "ii-E) 



+ — ^^-^^ (a2&2Mj + My(63g • Pj 



M IP 



/I 



27 



+aiihEfq- Pf) + a^iauEfMf + bsiMfq' + {q- PfY - a^q- PfEJ 
Ejq^ 



2M\Pj\ 



(3 cos^ ^ - 1) M;(a4 - 0264) +hEf{Mf + a^q- Pj - aua^Pf 



(C15) 



(fq Aaibi{M + Mf 



M] 



Ef q cos 9 



an) a2(&2^/ + hq- Pf) + Mj 



+a^hEf{q^ - aiiq- Pf)) + ^(cos' ^ " 1) {h{Mf + a^q- Pf - a^E))) 



(C16) 



d). For the charmonium he (J^^ = 1^~) and form factors the Vq, Vi, V2, V3 

The B-S wave function (solution of Eq. (IB9I) . Eq. ( JBIOI) and Eq. ( IBllI) under new method 
to solve the coupled equations) of h^ is read as: 



^h^Q') = (l'±-£ 



nw) + f2W)jf^ + fsW) ^± + f'M')^ 



M] 



75, 



with the constraint on the components of the wave function, 

f2iq')Mf 



/3(?0 = o, f'M) = 

Then we have the reduced wave function ip^^{q'): 
1 



m\ 



l + a^^ + aa^l 



{±11 
Mf 



75, 



a.^-^UlM^m')^) .a, „5 



03 = , 04 



w'l 



With Eq. ( llTl) . the form factors may be presented by overlapping integrations: 

d^q Aaibi q^ 



Vr 



fcos" ^ - 11 



(27r)3 M/(M + M/) 2 
aih2Ef + a2bzEf — hi^Mf + a^bi^q- Pf 



(C17) 



(Cl^ 



(C19) 



^1 



cP'q Aaihi 
(27r)3 M 



/ 



E 



f , q cos 9 , 






\Pf 



an 



0462-5'/ + aua^b^Efq- Pf + 62^/ + «2&3^- -P/j 



l^^hM' + {2E) + M^) cos^ 9) - «nSi^^ 
{aih2Ef + a2h-iEf - h^Mf + 0464^ ■ P/^ 



(C20) 



K. 



c/3g 4Mai6i 
(27r)3 Mf 



Ef qcos9 



) (—0464^^ + a2 — 0462-^/01] 



+ (aii ^^^T + ^^^-. (3 cos^ g - 1)) (04^2^/ + 02^3^/ - hMf + ttihq- P/)](,C21) 



28 



^'^~J (27r)3 MM ^|('°' ^ - ^) f^^(^2 + &4i?/«ii) + ha,] . (C22) 

e). For the charmonium Xc2 (J^^ = 2++) and form factors the k, ci, c,, h 

The B-S wave function (solution of Eq. ( 1B9I) . Eq. ( IBIOP and Eq. ( IBllI) under new method 
to solve the coupled equations) of Xc2 is read as: 



¥>xAg') = e,.qT{qT[niq') + ^J^iq') + j^Jsi^') + ^^^''^' 



M 



2 f'M')] 



+Y[Mfr,{q')+ rfaq')+ ^'Mq')] + -^fsiQ'V^'Pfaql^lsl,} , (C23) 

with the constraint on the components of the wave function: 



Mfm\ 



m'l 



where e^j, is a tensor for J = 2. Then we have the reduced wave function y^tj^i^ ^^' 



M 



(C24) 



/ 



V 



with 



fli 



gx ^^ I {f^{q')w'2-n{q')m',)Mf ^ ^^ U',{q')w', - f^{q')m',)Mj 
2Mf'm[ 2m[w2 ' 2m[w2 



1 n{q'')M] 



a3 = -ni + 

2 2m 1^2 



1 , w\ , 



2m'^w^ 



2 ^ rri'i ^ "^ ' 2m'iW2 

Mf Mfm\ M'l 

Q5 = ^^n2,a6 = ^\ n2, aj = , ag = —^n2 , 



2w[ 
With Eq. f fT7|) . the form factors may be presented by overlapping integrations: 



1 _. 777 1 _, _, in 

} in, } m. 



k 



d^q 



Ahi 



Q .__ 2 



cos^ ^ - 1) 



(27r) 3 M|(M + M/) 
-aiiEf{Mf{-a-i - 0363^/ + aih^Mj) + 04(62^/ + 64?- i^/)) 
-(aiias^s-El - a^hM'j + agfes? ■ -Pf) - aiiEf{Mf{-a^ - 0263-^/ + aibiMf) 
+ai{b2Ef + 645- Pj) + 2a8&3^/)) 
-E'/(aii T^T'' y^^^f + (^shEfq'^ + aQb2EfMf + af^h^Mfq- Pf - andshEfq- Pj 



\P 



Ell^S^^l _ cos^ e) ((Mj(-a3 - 0263^/ + a.hMf] 



+04(62^/ + &4g ■ i^/) + 2ashEf)) 



(C25) 
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Cl 



(27r)3 M| 



'^11 TT 1 <^iitt(~'^ii'^4^2-E/ " aua^b^Efq- P^ 



M 



M 



f 



+aib2MJ + a2hMfq- Pf + a^Mf{h^q^ + anEf - anh^q- Pf)) -2anashEfq- Pf/M 
E/qcos ( _^^^fY^Mf{-a^ - 0263^/ + aih^Mf) + 04(62^/ + h^q- Pf)) 



M\Pf 



M' 

Ef . . 

~aii — {anashEj - a^bsMj + aghq- Pf) 



E 



+2Q;ii-y (-0:110462 -E/ - Q;iia4fe4£'/g ■ Pf + 0162^^ + a2b^Mfq- Pf 
+a^Mf{hzq^ + auEf - auh^q ■ Pf)) - aiilash^Efq- Pf/M 
-an-^iauasbaEJ + a^h^^MJ + a^hq- Pf] 



+ 



q^ 



2M2|P^|2 



-Mj + {2E) + M}) cos^ 



M{a4b4q^ + a2Mf + 0:110462 -E/ — oiiOaM/ — og — a^h^Mf + onas&s-E/) 
-Q;iiE/(M/(-a3 - a2hEf + aih^Mf) + a4{b2Ef + 64^- P/)) 

-(aiifls^s-E/ - O'bhM'j + ag^ag ■ -P/) 
-Q;iiE/(M/(-a3 - a2hzEf + aih^Mf) + a4{b2Ef + 64^ ■ P/) 



+20863^/))) - 



g^i?/ cos ^ 

2M2|Pj|3 



(3M| - (2Ej + 3Mj) cos^ 



Mf{-a3 - a2hEf + a^b^Mf) + a4{b2Ef + b^q- Pf) + 208^3^/ 



(C26) 



C2 



(Pq AbiM 
(27r)3 Ml 

£^/g cos ^ 



an -J (^aiiE/(a464g^ + a2Mf + 0110462-^/ - Q;iia3M/) 



-a^b^q^ - aQb2Mf + aua^Ef) 



M\Pf 



~ai^Ef){Mf{-a^ - 0263^/ + aib^Mf) + a4{b2Ef + b^q- Pf)) 



-an{ana8hEJ - a^^bsMfasb^q- Pf) - anEf{a4b4q^ + a2Mf + 0:110462-^/ - aiia3M/) 

+ {-asb3q^ - aeb2Mf + auagEf) 
—aiiEf^a^b^q^ + a2Mf + 0:110462 -E/ — ana^Mf —ag — aeb^Mf + aua^b^Ef)) 

J2 



+ 



y / 71 7^2 I ( r\ rp2 I /i/f2\ „^^2 n\ ( ^ u ^2 



2M|PfP 



-Mf + {2Ef + Mi) cos^ 6) [a^biq^ + a2M/ + 0110462^, 



/I 



^/ 



-«iia3Mj — a^ — agb^Mf + 0110863-^/) 



2M|Pj|2 



(Scos^^-l) 



-aiiE/(M/(-a3 - 0263^/ + aibiMf) + 04(62^/ + b^q ■ Pf)) 
-{aiiagb^EJ - a^b^Mf + a8&3^- Pf) - anEf{Mf{-a^ - a2bsEf + aib4^Mf 
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q^ cos 9 

mPf\ 



+04(62^/ + b^q- Pf) + 2ashEf))) - ^^—^[{AEJ + Mj) cos^ 9 - {2EJ + Mj)] 



{Mf{-a^ - 0263^/ + aihMf) + 04(62^/ + hq ■ Pf) + 2a8&3^/)) ) ] , (C27) 

«ii) -jjiashq'^ + aQb2Mf + auasEf) 



, d^q AhAM + Mf 
h — ' 



(27r)3 M] 



, q cos 9 ^ ^ Ef 2 



\Pf\ ' M 



Ejq cos 9 Etq cos 9 o 2 

M\Pf\ M\Pf\ ^ ^ 

^ '-M) + {2E) + M)) cos^ 9)[as - aebiMf + auagbaEf) 



2MPf 

+ "^ "i 2 (3 cos^ - l){b3{aiiasE'i + aM - agq- Pf)) 
2MPf 



^.--.2. n f ^, Ef, 



-(cos 9-1) {-aii-j-{2{h'iMf{a^aii - 02) + 04(62 + aii&4-E/)) + 4a8&3 
+ ^^ ^T (1 - cos^ ^)(2(63M/(a3aii - 02) + 04(^2 + aub^Ef)) + 40363 



2M|P/|2 



. (C28) 
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